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' We perform a new study of the chances of the fluorescence detector (FD) at the 

. Pierre Auger Observatory to detect the tau leptons produced by Earth-skimming 

I ultra high energy z^t-'s. We present a new and more detailed evaluation of the effective 

' aperture of the FD that considers a reliable fiducial volume for the experimental set 

I ' up. In addition, we take into account the real elevation profile of the area near Auger. 

2 , We find a significant increase in the number of expected events with respect to the 

I predictions of a previous semi-analytical determination, and our results show the 

■ enhancement effect for neutrino detection from the presence of the near mountains. 
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Neutrinos constitute one of the components of the cosmic radiation in the ultra 
high energy (UHE) regime. Since we have detected ultra high energy cosmic 
rays (UHECR), the presence of a secondary UHE neutrino flux is guaranteed as 
a result of the vr-photoproduction, due to the interaction of hadronic UHECR 
with the cosmic microwave background. The detection of these cosmogenic 
neutrinos [1,2], in addition to a possible primary neutrino flux, would provide 
precious information on the physics and position of their powerful astrophys- 
ical sources. On the other hand, copious neutrino fluxes are also predicted 
in more exotic top-down scenarios where relic massive particles, produced at 
the first moments of the Universe, decay into UHE lighter particles, among 
which neutrinos and photons are expected. In any case, the detection of UHE 
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neutrinos would significantly contribute to unveiling the still unknown origin 
of UHECR. 

Due to the very low expected fluxes and the small neutrino-nucleon cross sec- 
tion, neutrinos with energies of the order of 10^^ eV and larger are hardly 
detectable even in the new generation of giant array detectors for cosmic radi- 
ation, like the Pierre Auger Observatory (Auger, in short) [3,4]. The detection 
of UHE neutrinos inducing inclined air showers was recently reviewed in [5]. 
In particular, a promising strategy concerning the detection of the tau leptons 
produced by Earth-skimming UHE i/^'s has been analyzed in a series of papers 
[6]-[19]. UHE r's, with energies in the range 10^^"^"^ eV, have a decay length 
not much larger than the corresponding interaction range. Thus, if a UHE 
crosses the Earth almost horizontally (Earth-skimming) and interacts in the 
rock, the produced r has a chance to emerge from the surface and decay in 
the atmosphere, producing a shower that in principle can be detected as an 
up-going or almost horizontal event. 

The aim of this letter is to perform a new, more refined, estimate of the effec- 
tive aperture of the fiuorescence detector (FD) at Auger to Earth-skimming 
UHE i/^-'s. A calculation of the number of possible up-going r showers de- 
tectable with the FD has already been performed in ref. [15] by using a semi- 
analytical computation. Our analysis represents a considerable improvement 
with respect to the estimate of this last work, since it uses a different method 
for calculating the number of u^/t events, which now includes a class of tracks 
neglected in the previous calculation. 

An additional improvement in our analysis comes from considering the effects 
of the topology around the Auger observatory site, by using a Digital Elevation 
Map (DEM) of the area around the experiment. A detailed DEM of the Earth 
surface is provided by ASTER (Advanced Spaceborne Thermal Emission and 
Reflection Radiometer) [20] which is an imaging instrument that is flying 
on Terra, a satellite launched in December 1999 as part of NASA's Earth 
Observing System. The available elevation map, GTOPO30, is a global digital 
model where the elevations are regularly spaced at 30-arc seconds. We show in 
Fig. 1 a 3D map of the relevant region around Auger. We will use this DEM to 
produce a realistic and statistically signiflcant sample of possible z^t/t tracks 
crossing the flducial volume of Auger, that will be used later to evaluate the 
real aperture of FD at Auger. 

We will define the Auger fiducial volume as that limited by the the six lateral 
surfaces (the subindex a — W, E, N, S, U and D labels each surface 
through its orientation: West, East, North, South, Up, and Down), and with 
Qa = {Oai 4>a) a gcucric direction of a track entering S^, as shown in Fig. 2. 
We have considered a simplified Auger area, given by a 50 x 60 km rectangle 
(an approximation to the real one, see ref. [4]), while the height of the fiducial 
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Fig. 1. A 3D map in longitude and latitude of the area around Auger with the 
elevation (not to scale) expressed in meters. The Auger position and surface, ap- 
proximated to a rectangle, is indicated in red. 
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Fig. 2. A simplified scheme of the Auger fiducial volume is represented (height not 
to scale). The lateral surfaces are labelled by their orientation. Two examples of 
entering tracks are also shown. 

volume was fixed to 10 km in order to be within the range of detection of the 
FD eyes. This is a conservative estimate, since we expect that the effective 
fiducial volume for the detection at the FD will be larger. 



Let be an isotropic flux of z/t- + i^t-. By generalizing the formalism developed 
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in ref. [13], the number of r Icptons emerging from the Earth surface with 
energy E^-, going through S„ and showering in the fiducial volume per unit of 
time (thus potentially detectable by the FD), is given by 



'""^^ -D fda. Ids.. fdE.. 



\ dt J ^ J J J dEjjdVta 

X j dEr eiEr) cos i0a) ka{E, ,Er;fa, ^a) , (l) 

where D is the duty cycle and e(£'^) is the detection efficiency of the FD, 
respectively [3]. The minimum energy for the r leptons, 10^^ eV, is chosen 
taking into account the energy threshold for the flourescence process [15]. The 
kernel ka{Ey jEr'jfa, ^a) is the probability that an incoming neutrino crossing 
the Earth with energy E,y and direction Qa, produces a lepton emerging with 
energy Er, which enters the fiducial volume through the lateral surface dSa 
at the position and decays inside this volume (see Fig. 2 for the angle 
definition). In Eq. (1), due to the very high energy of 1/^-, we can assume that 
in the process u-,- + N ^ r + X the charged lepton is produced along the 
neutrino direction. 

As already shown in details in ref. [15], this process can occur if the following 
conditions are fulfilled, 

1) the Ut- with energy E^ has to survive along a distance z through the Earth; 

2) the neutrino converts into a r in the interval z,z + dz; 

3) the created r emerges from the Earth before decaying with energy Et-; 

4) the T lepton enters the fiducial volume through the lateral surface at the 
point fa and decays inside this volume. 

1) The probabifity Pi that a neutrino with energy Ei, crossing the Earth 
survives up to a certain distance z inside the rock is 



where 



^CciEu) - ^yN(T^ \ ^ AT (3) 



is the charged current (CC) interaction length in rock [q^ ~ 2.65 g/cw?). 
A detailed discussion of an updated evaluation of the neutrino-nucleon cross 
section, aQ^{E^), can be found in ref. [15]. The dependence of X'^^j on the 
particular track direction can be safely neglected, since it is well known that 
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the interesting events are almost horizontal and thus the experienced Earth 
density is essentially^ equal to Qg- 



2) The probability for i/,- — > r conversion in the interval [z, z + dz] is 

^ , dz , , , 

^CC K^y) 



3) The probability Pc that a charged lepton survives as it loses its energy 
travelling through the Earth is described by the coupled differential equations 



'-r=-^P' . (5) 

dz C Tr tjr 

dE 

-^^-{Pr + lrEr) Qs ■ (6) 

Here = 1.77 GeV, ~ 3.4 x 10~^"^s denotes the r mean lifetime, whereas 
the parameters (3r — 0.71 x 10~^ cm^ g~^ and — 0.35 x 10"-^^ cm^ g~^ 
GeV"^, as discussed in ref. [15], fairly describe the r energy loss in matter (for 
further references and a recent discussion, see [21]). We denote the transferred 
energy as E^ = E^{Ey) = (1 — {ycc))Ey (see ref. [15] for details). By solving 
Eq.s (5) and (6) at the emerging point on the Earth surface one has 



^ _ I3tE^^{E^) CXp {-QsPr (^max - z)} 

" + 7, £;0(E,) (1 - exp {-Qs PriZr.^ ' z)}) ' ^ ^ 



where 



E°(E.)(/3. + 7.i^.) .._ ^rlr ... 



The quantity 2;max = -^max (^a , ^o) represents the total length in rock for a 
given track entering the lateral surface Eo of the fiducial volume at the point 
fa and with direction 

The energy Et of the exiting lepton must be consistent with Eq. (8). This 
condition is enforced by the presence of a (^-function in the final expression of 

^ We neglect the fact that UHE neutrinos may traverse a significant amount of water 
(as those from the West), which could enhance the arrival of UHE tau neutrinos at 
the largest energies, see e.g. [14]. 
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the probability P3, 



p ^ p x( p _ E^jEu) exp {-Qs /3r(^max - z)} \ 



4) Once the r lepton has emerged from the Earth surface, its showering prob- 
abihty is determined by the total distance it has to travel for reaching the 
fiducial volume. If we denote with Aout = Aout (^a , ^a) the total length trav- 
elled in the atmosphere by the r before reaching the point fa on the surface 
Efl, and with Ain = Ain {fa , ^la) the length of the intersection of the track with 
the fiducial volume, the decay probability inside the fiducial volume is given 

by 

n f Kutm-r\ ( f Ai„m^l\ 



Collecting together the different probabilities in Eq.s (2), (4), (10) and (11), 
we have 

•2 m ax 

ka{E, ,Er;fa,na)= J Pi P2 P3 A dz . (12) 



Since the flux can be reasonably assumed isotropic, namely 



dEy dVta dEj, dn 



Va , (13) 



we can rewrite Eq. (1), summing over all the surfaces, as 

where the effective aperture, A{Ei,) = Y.a'^aiEi,), is the sum of each surface 
contribution, 

Aa{E,) = J dErKa{E,,E,) , (15) 



and 



K,{E, ,Er)^l dQa I dSa COS {9a) e{Er) ka{E, ,E,;fa, Qa) 
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Fig. 3. The total effective aperture A{E,y) is plotted versus the neutrino energy 
(solid line). The dashed line corresponds to the same quantity as obtained in ref. 
[15] for H = 30 km. 

= J dfia J dSa cos {9a) e{E,) J P^P^P^PAdz . (16) 



In order to get the explicit expression for Ka{Ey , E^) an extremely involved 
integration has to be performed, which requires the computation of all the 
properties for each track (ra, Via), taking into account the DEM of the Auger 
site. To this aim a suitable approach is based on the following procedure: 
we use the available DEM of the Auger area to isotropically generate a large 
number of oriented tracks (let us say A^) which cross the Auger fiducial volume. 
If we denote with A^^ the subset of the N tracks which enter through surface 
E(j, then the kernel Ka{Ei, , Et) can be well approximated by the expression 

S 

Ka{E,,Er)^27re{Er)-^ E cos (^, J ka{E, , Er ■,n^ M ■ (17) 

We show in Fig. 3 the aperture A{E^) as a function of the neutrino energy 
(solid line), compared with the results of ref. [15] (dashed line). Remarkably, 
an enhancement factor approximately of one order of magnitude is found with 
respect to the semi-analytical results of ref. [15]. The origin of the difference 
between the two calculations can be explained in the following way. In ref. [15] 
an event is rejected if the r takes more than 30 km to decay; this means that 
almost horizontal events (where the r emerges from the Earth far from the 
lower surface of the experiment, travels more than 30 km before entering in the 
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Fig. 4. The effective apertures AaiE^) defined in Eq. (15) are plotted versus the 
neutrino energy. The thin sohd hne corresponds to the same quantity as obtained 
in ref. [15] for = 30 km. 

fiducial volume, and then decays inside it), preferred by the lateral surfaces 
due to the cos(^j^) term in Eq. (17), were not included in the calculation 
of ref. [15], while they are present in this analysis. In Fig. 4 the different 
contributions from each surface are plotted together with the results of ref. 
[15] (thin solid line). As one can notice, the contribution coming from the 
tracks going through the Down surface (thick solid line) is in nice agreement 
with our previous calculation. Moreover, it is possible to distinguish a North- 
South effect and, more consistent, a West- East effect. 

It is worth observing the different high energy behavior of the aperture cor- 
responding to the D surface, Aj^^ with respect to the others. The increase of 
the neutrino energy would select the almost horizontal tracks, which however 
are depressed, for the D-surface only with respect to the lateral ones, by the 
factor cos (6'^^ ) . 

We can use the expression in Eq. (14) to obtain the yearly number of r show- 
ering events at the FD of Auger (assuming a duty cycle D = 10%). In Table 
1 these rates arc reported for the same UHE neutrino fluxes considered in 
section 2 of ref. [15] (see in particular figs. 1 and 2), and described in a series 
of papers [22,23,24,25,26,27,28]. The three GZK fiuxes refer to three possible 
scenarios for cosmogenic neutrinos, which are those produced from an initial 
fiux of UHE protons. Instead the NH (New Hadrons) and TD (Topological 
Defects) cases are two examples of exotic models capable of generating the 
UHECR above 10^° eV, with large associated neutrino fiuxes. For each neu- 
trino fiux (fixed column), we list the total number of yearly events as well as 
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Table 1 

Yearly rate of Earth-skimming events at the FD for the different neutrino fluxes 
considered in ref. [15]. The number of r's showering into the fiducial volume that 
enter through each lateral surface are reported, as well as the total number of events 
for each flux. For comparison, we include the corresponding results from ref. [15]. 

the contributions from each lateral surface. 

Some comments are in turn. The number of events crossing the bottom surface 
is in fair agreement with the previous analytical result of ref. [15]. However, the 
total number of events is a factor 4-6 larger (depending on the model consid- 
ered), showing that the main contribution to the number of events is coming 
from almost horizontal showers, where the r emerges from Earth surface far 
away from the Auger fiducial volume and decays inside it. The enhancement 
factor depends on the different features of the fluxes used in the analysis. For 
example, for the three GZK models in Table 1, this factor ranges from 3.7 to 
6.2, corresponding to a hardening of the differential fluxes in energy (sec Figs. 
1 and 2 of ref. [15]). Instead, the enhancement is roughly the same for the last 
two models in Table 1 despite the fact that they have a different spectrum in 
the high energy range, due to the suppression of the very high energy neutrino 
events which escape without showering. 

As one can see from Table 1, a significant difference in the number of events ex- 
ists between the Surfaces W (facing the Andes) and E, which shows a mountain 
effect. This enhancement is mainly due to the largest amount of rock encoun- 
tered by horizontal tracks coming from the west side. If we define, as a measure 
of the effect of mountains, the difference of the number of expected events en- 
tering the volume from W-surface and E-surface, divided by their average, the 
effect can be quite remarkable (even order 30%). This effect would also be 
larger by comparing the exclusive apertures, but in this case the difference in 
the apertures, which is larger for larger neutrino energy, should compel the 
fast decreasing with energy of neutrino fluxes. Of course, on the total number 
of events this 30% effect is diluted because it concerns only one surface among 
four lateral ones. Moreover, a similar but smaller effect (of order ~ 20-25%) is 
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0.043 


0.116 
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Table 2 

Yearly expected number of down-going events at the FD, due to the showering of 
Vf. and Vt inside the Auger fiducial volume. The different predictions refer to the 
same fluxes of Table 1 and to a zenith angle larger than 60° . 

also present in the difference between the number of events from the Surfaces 
N (facing the higher part of the cordillera) and S. 

In conclusion, the largest contribution through the surfaces E, W, S and N 
makes the total number of yearly Earth-skimming events larger than the previ- 
ous estimates of ref. [15]. This in turn increases the detection chances of UHE 

Ur at the FD of Auger, which seem realistic even for conservative neutrino 
fluxes like GZK-WB, considering that data will be taken over many years of 
observation. 

In analogy with the analysis of ref. [15], we report in Tabic 2 the expected 
numbers of yearly down-going events at FD, produced by UHE and 
showcring inside the fiducial volume ^ . The numbers arc calculated in a similar 
way to the Earth-skimming case, but with a different kernel in Eq. (12), 

ka{E,^ ,E^;f,,n,)^5 (e^ - El{E,S) Pi'' P2" (18) 



where a = e,T and E^{E^^) is the transferred energy to the charged lepton. 
For a i^e we have 



which take into account the probability that the neutrino survives until the 
fiducial volume producing an electron which initiates a shower inside the de- 
tector. Instead, for a Vr we have 



^ From the study of t'ct'^-induccd shower in the atmosphere of ref. [29] it is argued 
that high energy f^'s have a very small probability of inducing a shower before 
reaching the ground. 
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Pi''^ = 1 - exp 



CT-r Er 



dz 



exp {-z/X-,,,iE,,:)} exp {-(A,„,, - .: ) in,/{rT, 



exp {-Aout/A^g(-g^)} - exp {-Aput mr/jcTr E^)} 



(20) 



which take into account the probabihty that the neutrino produces a r which 
survives until the fiducial volume and initiates a shower inside the detector. 



In this letter a new and more detailed computation of the aperture for Earth- 
skimming Vr of the fluorescence detector at the Pierre Auger Observatory 
has been performed. The evaluation has been carried out by using a different 
approach with respect to the previous semi-analytical one of ref. [15] and 
taking into account the real elevation profile of the area around Auger. The 
obtained results show an increase of the effective aperture, and correspondingly 
of the number of the expected events. This larger result is mainly due to 
the contribution of almost horizontal Earth-skimming r's, emerging from the 
Earth surface far away from Auger but decaying inside its fiducial volume. 
Remarkably the presence of mountains near Auger leads to an enhancement 
factor on the total number of events of the order 1.3-1.6. As previously shown 
(see e.g. the discussion in [5]), our results indicate that the number of UHE 
neutrino events at Auger is comparable for showers induced by down-going 
UHE neutrinos and by the decays of r's produced by Earth-skimming UHE 

l^T-'S. 



It is worth observing that the efficiency of the FD detector, whose parame- 
terization we have used in Eq. (1), can be considered as oversimplified and 
that one should expect a behaviour of FD much more complicated with en- 
ergy and/or geometry of the event. Unfortunately, a complete and exhaustive 
analysis of this issue is still absent and thus not available for any analysis like 
ours. Nevertheless, our predictions can be considered as a first but probably 
useful step in the right direction, which should motivate more study inside the 
Auger collaboration for example in order to obtain a reliable efficiency of FD 
with energy and for horizonal showers. 



Possible future improvements of the present work should take properly into 
account the response of the flourescence detector, as well as the real geometry 
of the Auger area and the fact that some of the UHE u^-/ r tracks could partly 
traverse water. 
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